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We study the non-equilibrium dynamics of a coherently split one-dimensional (Id) Bose gas by 
measuring the full probability distribution functions of matter- wave interference. Observing the 
system on different length scales allows us to probe the dynamics of excitations on different energy 
scales, revealing two distinct length-scale dependent regimes of relaxation. We measure the crossover 
length-scale separating these two regimes and identify it with the prethermalized phase-correlation 
length of the system. Our approach enables a direct observation of the multimode dynamics char- 
acterizing one-dimensional quantum systems. 



The non-equilibrium dynamics of many-body quan- 
tum systems and their pathway towards equilibrium is 
of fundamental importance in vastly different fields of 
physics. Open questions appear, for example, in high- 
energy physics for understanding quark-gluon plasma [1- 
3], in cosmology for describing preheating of the early- 
universe [4], or in the comprehension of relaxation pro- 
cesses in condensed matter systems [5, 6]. 

Due to their isolation from the environment and their 
tunability, ultracold atom systems have triggered many 
studies of non-equilibrium dynamics in closed interact- 
ing quantum systems, with a particular interest drawn 
towards quantum quenches [7, 8]. Important questions 
are related to systems where the dynamics is constrained 
by several constants of motion [9], and to the possible 
description of non-equilibrium states by generalized sta- 
tistical mechanics ensembles [10, 11]. 

Recently, we reported the experimental observation of 
prethermalization in a coherently split one-dimensional 
(Id) ultracold Bose gas [12], made possible by a charac- 
terization of the dynamical states through measurements 
of full distribution functions [13, 14]. Prethermalization 
[15] was understood as the rapid relaxation to a steady 
state exhibiting thermal-like properties but differing from 
the true thermal equilibrium that is eventually expected 
to occur on longer time-scales [16-20]. 

In this letter, we study the relaxation process [21] lead- 
ing to the prethermalized state by measuring the full 
distribution functions (FDFs) of phase and contrast of 
matter-wave interference. We probe the Id system on 
different length scales to investigate its multimode dy- 
namics which reveals two distinct regimes separated by 
a characteristic crossover length-scale. We measure this 
characteristic length-scale and identify it with the effec- 
tive thermal phase-correlation length of the prethermal- 
ized system. 

We prepare a quasi- Id Bose gas of several thousand 
87 Rb atoms in an elongated magnetic microtrap on an 
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FIG. 1: (Color online) Matter- wave interferometry and cor- 
relation measurements, (a) The experiment is initiated by 
coherently splitting a single quasi- Id Bose gas into two un- 
coupled gases using a horizontal double- well potential. The 
system is then held to evolve for a variable evolution time t e , 
before being released from the trapping potential. The con- 
trast of the resultant interference pattern is a direct measure 
of the relative phase fluctuations (b). Integrating the inter- 
ference pattern over a length L and fitting the resulting line 
profile (c) gives the phase &(L) and the contrast C(L) which 
are plotted as circular statistics (d). (e) The whole process 
is repeated many times (typically 150) to map the FDF for 
a particular length scale L and evolution time t e . The point 
plot in (e) is then converted into a density plot (f). 



atom chip [22] at a (tunable) temperature between 20 
and 120 nK. The gas is coherently split along the ra- 
dial direction using a symmetric radio-frequency dressed- 
state double- well potential [23], creating two uncoupled 
Id gases (Fig. 1(a)). The system is then held in the 
double-well trap for a variable evolution time t e before 
being released and allowed to fall under gravity during a 
time-of- flight of 16 ms. The resulting matter- wave inter- 
ference pattern (Fig. 1(b)) is recorded using absorption 
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FIG. 2: (Color online) Multimode dynamics revealed by FDFs of matter-wave interference. The probability density of contrast 
C(L) and phase &(L) of interference patterns is measured for each integration length L (horizontal axis) and evolution time 
t e (vertical axis). Red and blue denote high and low probability respectively and the color-map is rescaled for each plot. 
For each value of L, the right (left) columns correspond to the experimental data (theoretical calculations), and full cloud to 
L = 100 /im. At t e = 1.5 ms the high contrasts and small phase spreads demonstrate the coherence of the splitting process. As 
time evolves, a steady state emerges and two distinct length-scale dependent regimes appear: the phase diffusion regime (A) 
and the contrast decay regime (B). For short (long) L, the phase is random and the probability of observing a high contrast 
is high (low), resulting in a ring (disk) shape in the density plot. The theoretical calculations take into account the technical 
noise on the relative atom number between the two wells (standard deviation measured to be 5% of the total atom number) 
and the error associated with the fitting of the interference patterns. 



imaging [24] and contains information about the local rel- 
ative phase (j)(z) = 0i(z) — 02 (z) between the two halves 
of the system, with Oi^(z) being the fluctuating phase 
profile of each individual gas. 

In Id systems, spatial fluctuations arising from excita- 
tions at different wavelengths strongly affect the physics 
[25], and probing the system on variable lengths L rep- 
resents a filter for the effects of these excitations [26]. 
In our experiment, we integrate the interference pattern 
longitudinally over a length L and extract a line pro- 
file, from which a contrast C(L) and a phase $>(L) are 
obtained (Fig. 1(c)). Repeated realizations allow us to 
measure the full probability distribution function (FDF) 
of phase and contrast for different evolution times t e and 
integration lengths L. Our measured FDFs are presented 
in Fig. 2. 

For the shortest evolution time (1.5 ms in Fig. 2), 
we observe high contrasts and small phase spreads for 
all lengths L, which demonstrates the coherence of our 
splitting process. After splitting, the phase fluctuations 
corresponding to the different excitations grow in magni- 
tude, resulting in a scrambling of the relative phase field 
4>(z). These spatial fluctuations of (j)(z) manifest itself 



in a randomization of the phase in the interference pat- 
terns, and in a decrease of the contrast which depends 
on the probed length-scale L. For short L, the sparsely 
populated high momentum modes satisfying k > 2tt/L 
do not lead to contrast reduction, resulting in a ring-like 
shape of the FDFs; this is the phase-diffusion regime (A 
in Fig. 2). For long L, many modes satisfying k > 2tt/L 
are populated and their dynamics leads to a scrambling 
of (j)(z) within the probed integration length, resulting in 
a significant reduction of the probability of observing a 
high contrast and in a disk-like shape of the FDFs; this 
is the contrast-decay regime (B). 

At longer evolution times (t e > 10 ms), the 
FDFs do not change significantly any more, reveal- 
ing the emergence of a steady state characterized by a 
crossover length-scale between the phase-diffusion (A) 
and contrast-decay (B) regimes. In our previous work 
[12], we identified this steady state as a prethermalized 
state. 

We now quantitatively analyse our observations in the 
framework of an integrable Luttinger Liquid (LL) theory. 
To describe the dynamics, we consider the LL Hamilto- 
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the time evolution of the operators </>(z) and h(z) rep- 
resenting the relative phase and the relative density of 
the system, respectively [21, 26]. Here, c = ^/pgJm is 
the sound velocity, K = n^p the LL parameter charac- 
terising the strength of the interactions and £ = Ti/mc 
the healing length; p is the Id density in each half of 
the system, g = 2huj±a s the Id coupling constant, m the 
mass of the 87 Rb atom, a s the scattering length and u±_ 
the radial trapping frequency [27]. The time evolution 
of n(z) and <j>(z) can be described in Fourier space by a 
set of uncoupled harmonic oscillators of collective modes 
with momentum k and energies ftuok ~ hc\k\, where the 
LL Hamiltonian is diagonal [14]. 

The rapid splitting process prepares a coherent super- 
position of the atoms in the two wells, resulting in high 
relative density fluctuations (n|)t=o = p/2 and small rel- 
ative phase fluctuations ((j) 2 )t=o = 1/2/9 [14, 27]. In the 
experiment, we observe enhanced global phase fluctua- 
tions for t e = 1.5 ms with respect to quantum noise, 
which are attributed to technical noise on the relative 
atom number between the two wells. For the dynamics 
of the contrast, the contribution of this technical noise is 
however negligible compared to that of the initial density 
fluctuations, for the probed integration lengths {L ^> £). 
Focusing on the contrast thus allows us to directly ob- 
serve the effect of quantum noise associated with the 
splitting process through the evolution of the many-body 
system. 

The results of the LL calculations, taking into account 
our technical noise, are presented in the left columns of 
Fig. 2. For all length scales and evolution times, we find 
very good agreement between the experimental and the- 
oretical FDFs. Our approach gives a direct and intuitive 
visualization of the multimode nature of the dynamics in 
Id Bose gases [28]. 

To study the relaxation of the contrast in more de- 
tail, we compare in Fig. 3 the measured mean values of 
the integrated squared contrast (C 2 (L,t e )} to the pre- 
diction of the LL theory, taking into account the imper- 
fections of our imaging system [27]. We observe that 
the characteristic time for the decay of the contrast is 
longer for long L than for short L, which is well cap- 
tured by our model. For a more quantitative compari- 
son, we fit the experimental and theory data by an ex- 
ponential decay. For the experiment, we find charac- 
teristic decay times r exp (L = 6 pm) = 0.53 ± 0.33 ms 
and r exp (L = 100 pm) = 2.89 ± 0.36 ms (standard error 
on fitted parameter). Fitting the theoretical calculations 
for our measured LL parameter K = 58 ± 4, we find 
r th (L 6 pm) 0.40 ± 0.11 ms and r th (L 100 pm) 
3.15 ±0.21 ms, in good agreement with the observations. 
The dependence of the relaxation time-scale on L is an 
additional clear signature of the multimode nature of the 
dynamics [14]. 
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FIG. 3: (Color online). Contrast dynamics. Measured values 
of the mean squared contrast for various integration lengths, 
corrected for the contrast reduction factor due to the imper- 
fections of our imaging system [27]. From top to bottom: L = 
18, 40, 60, 100 pm. The lines show the results of the Luttinger 
Liquid calculations for these integration lengths. The theory 
data for all L have been rescaled by the same factor r = 0.74 
[29]. Apart from this factor, no fitting parameter is used. We 
observe a relaxation process in which a steady-state is estab- 
lished on a time-scale depending on L and corresponding to 
the dephasing of the different excitations probed within that 
L. 



Fig. 3 also reveals the onset of the prethermalized state 
for t e > 10 ms which is characterized by high (low) con- 
trasts at short (long) integration lengths, in line with the 
probability density plots A (B) of Fig. 2. In the rest of 
this letter, we investigate the crossover between the high 
contrasts and low contrasts regimes observed in Figs. 2-3. 

Within the LL theory, the crossover length-scale sep- 
arating the phase diffusion and contrast decay regimes 
can be calculated analytically [14]: Iq = 2% 2 /mg = 
h/ma s uj±. For the data in Figs. 2-3 we measure uj±/2tt = 
1.40 ± 0.08 kHz and calculate l = 15.8 ± 0.9 pm. To 
understand the nature of this crossover length-scale, we 
consider the decay of the mean squared contrast (C 2 (L)) 
as a function of L for t e > 10 ms (inset of Fig. 4). To fit 
this decay, we use our knowledge of the correlation func- 
tions in the prethermalized state, which have been shown 
to be thermal-like [12, 14]. Using the direct link between 
the mean squared contrast and the relative-phase cor- 
relation function [27, 30], a fit of (C 2 (L)} allows us to 
estimate the effective thermal phase-correlation length 
A e ff of the system. For the data in Fig. 4 (inset) we find 
A e ff = 16.6 ± 0.9 pm, close to the value of Iq calculated 
above using the independently determined experimental 
parameters. Note that our method of estimating the cor- 
relation length does not require precise measurements of 
the system parameters to fit the FDFs, but only needs 
the knowledge that the correlation functions are effec- 
tively thermal. 

The equivalence of the crossover length-scale Iq and 
the effective thermal phase-correlation length A e ff can be 
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FIG. 4: (Color online) Comparison of the relative-phase cor- 
relation length of the prethermalized system (A e ff) and of a 
system at thermal equilibrium (A^). Blue squares: values of 
A e ff obtained by varying the temperature Tinit before split- 
ting the quasi-condensate (see main text). Error bars denote 
one standard deviation and are obtained by a bootstrapping 
method [27]. The red shaded area corresponds to the 95% 
confidence interval around the theoretical calculation (line) 
accounting for the uncertainty on the experimental param- 
eters. Green circles: correlation length A^ of a system of 
two quasi-condensates in thermal equilibrium, normalized to 
the mean of the different densities. Gray shaded area: 95% 
confidence interval around the theory (dashed line). Inset: 
estimation of the correlation length. The rilled circles are the 
values of (C 2 (L)) in the prethermalized state (data of Fig. 3 
for t e > 10 ms), normalized to the value at L = 6 pm. The 
solid line is a fit to an effective thermal equilibrium theory 
with the correlation length being the only free parameter. 

understood in the following way. In the case of two un- 
coupled quasi-condensates in thermal equilibrium with 
temperature T, the thermal phase-correlation length is 
given by A^ = h 2 p/mk&T [31, 32]. Due to our rapid 
splitting process, the energy initially stored in the sys- 
tem is equally distributed between the different /c-modes 
of the system [12, 14], resulting in thermal-like correla- 
tions characterized by an effective temperature k^T e R = 
(H)\t=o = pg/2- The effective thermal phase-correlation 
length can thus be identified with A e ff = Ti 2 p/mk^T e R = 
2ft 2 /mg, and is equivalent to the crossover length scale 
Jo- 

While the prethermalized system reveals thermal-like 
correlations, its correlation length A e ff depends only on 
the Id coupling constant g, in contrast to a system of 
two quasi-condensates at thermal equilibrium where A^ 
is a function of density and temperature. To reveal this 
difference experimentally, we varied the initial temper- 
ature Ti n it of the quasi-condensate before splitting and 
measured the correlation length by fitting the decay of 
(C 2 (L)} (see inset of Fig. 4). The temperature before 
splitting, Ti n it, was obtained through measurements of 
the second order correlation function of longitudinal den- 



sity fluctuations after time-of- flight [27, 33]. The results 
(blue squares in Fig. 4) show the independence of A e ff 
from Tinit- When varying T^it, we further observed that 
the evolution of the FDFs remained close to that pre- 
sented in Fig. 2. In particular, the steady states reached 
in the evolution exhibited the same crossover between 
the phase-diffusion and contrast-decay regimes. These 
observations confirm our interpretation of the crossover 
length-scale Iq as the prethermalized phase-correlation 
length A e ff of the system. Performing experiments with 
two quasi-condensates in thermal equilibrium with tem- 
perature T and prepared in the same double-well trap 
[27], we observe A^ ex 1/T (green circles in Fig. 4), em- 
phasizing the different characteristic scalings of the ther- 
mal and prethermalized phase-correlation lengths. 

In conclusion, we demonstrated the multimode nature 
of the dynamics of a coherently split Id Bose gas using 
matter- wave interferometry. Our approach allows the di- 
rect observation of the effect of the quantum noise present 
in the initial non-equilibrium state on the evolution of 
the many-body system. We showed the emergence of 
two distinct regimes of relaxation separated by a char- 
acteristic crossover length scale, which we identified as 
the prethermalized phase-correlation length of the sys- 
tem. This characteristic length-scale, 2ft 2 /mg, reflects 
the parametrization of the many-body theory describ- 
ing the dynamics of the system by a single parameter, 
g. For our closed system to thermalize and the phase- 
correlation length to reach the true thermal equilibrium 
value A^, the integrability of the LL theory needs to be 
broken. Investigations into the processes leading to the 
full thermalization of the system are ongoing. 
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